> 

(N 



in 
o 



X 



Acta Phys. Hung. A 19/1 (2004) 000-000 

' HEAVY ION 

PHYSICS 



Numerical Simulation of Non-Abelian Particle-FTeI3~ 

^ ; Dynamics 
O . 

CN . Adrian Dumitru and Yasushi Nara 
> , 

Q _ Institut fiir Thcorctischc Physik, Johann Wolfgang Goethe Universitat, 

' Max von Laue Str. 1, D-60438 Frankfurt am Main, Germany 

o , 

^ Received 21 November 2005 



Abstract. Numerical 1D-3V solutions of the Wong- Yang-Mills equations with 
anisotropic particle momentum distributions are presented. They confirm the 
existence of plasma instabilities for weak initial fields and of their saturation 
■ at a level where the particle motion is affected, similar to Abelian plasmas. 

The isotropization of the particle momenta by strong random fields is shown 
explicitly, as well as their nearly exponential distribution up to a typical hard 
scale, which arises from scattering off field fluctuations. By variation of the 
lattice spacing we show that the effects described here are independent of the 
' UV field modes near the end of the BriouUin zone. 
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1. Introduction 

Recently, it has been realized that non- Abelian collective plasma processes such as 
Weibel-like instabilities might play an important role for the thermalization process 
in the early stage of high-energy heavy-ion collisions. This was the central topic 
of the workshop on "Quark-Gluon-Plasma Thermalization" in Vienna [1]. If so, a 
quantitative understanding of such processes will be crucial to answer, for example, 
the question about the maximum temperature achieved in such collisions at the 
BNL-RHIC and CERN-LHC colliders. 

The physics of non- Abelian plasma instabilities in the context of rclativistic 
heavy-ion collisions has been discussed in some detail in a recent review [2] and in 
many contributions to this workshop. We shall therefore refrain from a detailed 
presentation here. Rather, we focus on illustrating the generalization of Abelian 
particle-in-cell simulations to the SU(2) gauge group. These provide some addi- 
tional insight into the physics of non-Abelian plasmas beyond the "Hard Loop" 
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approximation which underlies much of the present analytical and numerical under- 
standing of SU(2) instabilities [3]. Simulations of the non-linear Vlasov-Yang-Mills 
theory account for the back-reaction of the fields on the particles, which damps (and 
eventually shuts off) the exponential growth of the chromo-magnetic fields. They 
also enable us to actually look at the time-evolution and eventual isotropization of 
the particle momenta themselves. Another motivation for performing full particle- 
field simulations is the potential interest in initial conditions where the fields are 
strong and immediately affect the particle motion. Some of our results have been 
published in ref. [4] . 



2. Particle-in-cell simulations for non-Abelian gauge theories 

We consider the classical transport equation for hard gluons with non-Abelian color 
charge in the coUisionless approximation [5,6], 

P^[d^ - gQ^'F^.d; - gfabcAlQ-dQ.]f{x,p, Q)^0. (1) 

where / denotes the one-particle phase-space distribution function. We employ the 
test-particle method, replacing the continuous distribution f(t,x,p,Q) by a large 
number of test particles: 

fit, x,p, Q) = -J— V 5{x - x,{t)){2nfS{p - p,{t))5{Q - Q,{t)) . (2) 

fi {t) , Pi {t) , Qi (t) are the phase-space coordinates of an individual test-particle. (We 
consider particles in the adjoint representation of color-SU(iVc), hence Q is a vector 
in — 1 dimensional color space.) This Ansatz leads to Wong's equations [5,6] 

^ = v., ^=gQnE''+v,xBn, ^=^gv^[A,,Q,,l (3) 

for the z-th (test) particle, coupled to the Yang-Mills equation 

-^ = E\ — = i?,F^"^-^^gfet.M(a;-a;fc), {^ = x,y,z), (4) 

in the temporal gauge = 0. This set of equations reproduces the "hard thermal 
loop" effective theory [6] near equilibrium. In the following, we assume that the 
fields only depend on time and on one spatial coordinate, x, which reduces the 
Yang-Mills equations to l-fl dimensions. The particles are allowed to propagate in 
three spatial dimensions. This is referred to as 1D-I-3V simulations. 

Numerical techniques to solve the classical field equations coupled to colored 
point-particles have been developed in Ref. [7]. Our update algorithm is closely 
related to the one explained there, which we briefly summarize. We employ the 
so-called Nearest-Grid-Point (NGP) method which simply counts the number of 
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particles N{j) within a distance ±a/2 of the jth lattice site to obtain the density 
rij = N{j)/a (with a the lattice spacing). If a particle crosses a cell, a current 
is generated. For example, if a particle crosses from site i to « + 1, 

Mt,^) = ^^S{'--'-^). (5) 

The color charge then has to be parallel transported to the next site, 

Q(* + l) = C/t(z)gW[/,(*) . (6) 

The gauge links arc related to the continuum fields via Ux{i) — cxp(igaAx{i)). 
In this way, the continuity equation for color charge is satisfied locally, together 
with Gauss's law. At icross we also update the particle momentum p^. by imposing 
energy-momentum conservation in the presence of the chromo-electric field E^ii). 
On the other hand, the rotation of a particle's momentum due to the color-magnetic 
field is updated in every time step. 

For 1D-I-3V simulations a major simplification arises from the fact that the 
transverse current 

J^{t,x) = -^y^Q,v^5{x-x,{t)) , (7) 

I 

can be updated continuously in time. Note that the color rotation due to the gauge 
fields Ay and (which in fact become adjoint scalars in ID) is also continuous in 
time. Therefore, our transverse current is very smooth and much less noisy than 
which is obtained in the impulse approximation. However, for 1D-3V simulations 
the longitudinal current can also be made sufficiently smooth by employing a large 
number of test-particles. Such a "brute-force" approach is no longer feasible for 
3D-3V simulations. 

To check the numerical accuracy, we have first performed simulations for isotropic 
momentum distributions, varying the lattice spacing and the number of test-particles 
(Fig. 1). The field energy density is determined from the lattice field strength 
El s 9o?E as g^e ~ (l/2)£j£/a^, plus the magnetic contribution. One observes 
that the time evolution stabilizes with increasing number of test-particles and de- 
creasing lattice spacing. However, to reach sufficient accuracy our runs required 
several hundred to a thousand test-particles per lattice site, corresponding to w 6 
hours run-time on a single-processor 2.4 GHz Opteron workstation per initial con- 
dition. It is therefore clear that a simple-minded extension of the point-particle 
algorithm to 3D-3V simulations is impossible. For multi-dimensional simulations, 
a generalization of current smearing from U(l) [8] to non-Abelian gauge groups is 
essential. Non-Abelian simulations which include fluctuations of the fields in the 
transverse plane would be very interesting because of indications that this leads to 
the development of a turbulent cascade which transfers energy from the soft un- 
stable modes to stable UV field modes (near k ~ 1/a). This process, which is due 
to the self-interaction of the gauge field in the SU(2) theory, effectively tames the 
exponential growth of the fields [9]. 
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Fig. 1. Time evolution of the average field energy density for isotropic particle 
momentum distributions on two different lattices (with the same physical size L) 
and for varying number of test particles. Nc = 2 color simulations. 



3. Anisotropic initial distribution 

In what follows, wc consider anisotropic initial momentum distributions of the hard 
gluons, 

f{p, x) oc cxpi-^pl + pl/p\,i,,d)5{px) . (8) 
This represents a quasi-thermal distribution in two dimensions, with "temperature" 

= Phard- 

The initial field amplitudes are sampled from a Gaussian distribution with a 
width tuned to a given initial energy density. We solve the Yang-Mills equations in 
A° = gauge and also set A = at time i = 0; the initial electric field is taken to 
be polarized in a random direction transverse to the x-axis. This initial condition is 
convenient because Gauss's law DiE'^ = p then implies local (color) charge neutrality 
at t = 0, allowing for a straightforward initialization. Of course, magnetic and 
longitudinal electric field components quickly build up as time progresses. 

We first show results for a relatively large separation of initial particle and field 
energy densities which should qualitatively resemble the conditions studied in [3,9]. 
The results shown in Fig. 2 correspond to a lattice of physical size L = 40 fm 
and = 1024 sites (the plots shown in ref. [4] correspond to the same L but 
half the number of lattice sites). The hard scale was chosen as phard — 10 GeV 
(in lattice units, pL,hard = aphard), and the particle density g^n = 10 fm""^ (we 
define the density in lattice units as tll = g^a^n). The above definitions of the 
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Fig. 2. Time evolution of the kinetic (particle), and magnetic and electric field 
energy densities for U{1) and 5(7(2) gauge group, respectively. 



lattice hamiltonian, fields and phase-space distribution function remove any explicit 
reference to the gauge coupling g from the lattice theory. 

For the Abelian theory we observe a rapid exponential growth of the magnetic 
field energy density, starting at about t/L p» 0.1; we repeat that in order to avoid 
a "fake" growth of the fields during this initial transient time, one has to ensure 
that the number of test particles is sufficiently large. At a time t/L « 0.4 the 
magnetic field strength has grown by about one order of magnitude. The fields 
clearly affect the particles, which loose energy. In turn, at this time the exponential 
growth of the magnetic fields is slowed down. The electric field grows less rapidly 
and equipartitioning is not achieved within the depicted time interval. 

The non- Abelian case features a rather similar evolution for short times {t/L/Nc ~ 
0.2 for electric fields and ;=! 0.3 for magnetic fields, respectively). We scaled time by 
1 /Nc because such a scaling is natural in the linear regime [3] . The growth of the 
magnetic field then saturates somewhat earlier than for the U(l) theory, and due to 
commutators, the electric field has more strength by the end of the simulation. Due 
to the somewhat earlier saturation of the instability, the colored particles loose less 
of their energy to the fields than was the case for electric charges. Nevertheless, at a 
purely qualitative level the U{1) and SU{2) simulations are not extremely different, 
which is due to the phenomenon of "Abelianization" in 1D-3V simulations [2,3]. 
This does not occur in 3D. 

In Fig. 3 we compare results obtained on two different lattices with = 512 
and Nx — 1024 sites, respectively, for the same set of physical parameters. One 
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Fig. 3. Time evolution of the field energy density on two diff'erent lattices (with 
the same physical size L); Np denotes the number of test particles per lattice site. 

observes that the growth rate of the instability, the saturation level and time are 
nearly the same. This confirms the underlying physical picture that the dynamics 
is dominated by the unstable soft field modes rather than UV-modes near the end 
of the BriouUin zone. If field modes with k ^ \/a would affect the dynamics then 
the continuum limit would not exist. 

Our 1D-3V simulations thus clearly confirm the existence of instabilities in non- 
Abelian plasmas and the idea of "Abelianization" , namely that the field growth is 
perhaps damped by self-interactions but does not shut down until the fields have 
grown so much as to affect the motion of the particles. Nevertheless, the number 
of e-foldings by which the field energy density grows is much less spectacular in our 
simulations than for simulations within the "hard-loop" approximation [3]. This is 
due to the fact that our initial field amplitudes are already relatively large (non- 
linear regime). At a technical level, point-particle simulations are not very well 
suited to study the extreme weak-field regime, which would require a prohibitively 
large number of particles. 

Once the fields have grown strong, they deflect the particles from their straight- 
line trajectories and finally lead to isotropization of their momentum distributions, 
shown explicitly in rcfs. [4]. This process represents the dominant contribution to 
the build-up of longitudinal pressure. Figs. 4,5 depict the evolution of the particle 
distribution function. This result was obtained with "strong-field" initial conditions: 
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Fig. 4. Initial and final particle dis- 
tribution functions for the strong field 
case on a ~ 512 lattice. 



10 

10=^^ 
10 

1 

10-^ 
10-=^ 
10"^^ 
10' 



g2f(p ) jnjtlal 
g'^f(p ) initial 

g^Kp'l final 
g2f(pp final 

thermal 




N,=1024, N =1000 

„j , I \ , , i__ 

1 



p /p^ , or p /p^ ^ 

■^x "^hard y "^hard 
Fig. 5. Same as Fig. 4 for a iV^ = 1024 
lattice. 



Phard = 1 GeV and initial field energy density « 10~^ GeV/fm'^/g^; the time scale 
is set by the lattice size, L = 10 fm for this simulation. When the separation 
between hard and soft modes is not so large, strong instabilities can not develop as 
the system approaches isotropy very quickly [4]. 

In fact, propagation in strong random fields not only leads to isotropic but even 
to nearly exponential particle momentum distributions, as can be seen from Fig. 4. 
All particles with momenta up to ~ Phard appear to be more or less thermalized. 
Once again, we check the dependence on the lattice spacing by comparing results 
obtained on different lattices. We confirm numerically that the process does not 
appear to be dominated by the ultraviolet modes of the fields on the lattice. 

Field fluctuations generate an effective collision term mediating soft exchanges. 
Defining 

}{x,p,Q) = {,f) + 5f, Al = {Al) + 5Al, (9) 

where ( ) denotes the ensemble average, and {5f) — {SA^j) = 0, one can obtain the 
Balescu-Lenard collision term from the fluctuation part, showing the correspondence 
between fluctuations and collisions in an Abelian plasma [10]. For the non-Abelian 
case, see refs. [11]. We also note that the mean entropy density is no longer conserved 
in the presence of fluctuations. 
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